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High energy amplitude in the dipole approach 
with Pomeron loops: asymptotic solution 



gested. This solution has several unexpected features: (i) the asymptotic amplitude is a 
function of dipole sizes and, therefore, this amplitude shows the gray disc structure at high 
energy, instead of black disc behaviour which was expected; (ii) the amplitude approaches 
the asymptotic limit in the same way as the solution to the Balitsky-Kovchegov equation 
does (oc exp(— CY 2 ) ), but the coefficient C in eight times smaller than for the Balitsky- 



high energy asymptotic behaviour by changing the initial condition from Z(Y; [v,i = 1]) = 1 
to Z(Y; [ui = 1 — 7o,i]) = 1- The value of 70 is determined by the process of merging of two 
dipoles into one. With this new initial condition the Balitsky-JIMWLK approach describes 
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In this paper an analytical solution for the high energy scattering amplitude is sug- 
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Kovchegov equation; (iii) the process of merging of two dipoles into one, only influences the 




the high energy asymptotic behaviour of the scattering amplitude without any modifications 
recently suggested. 
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1 Introduction 



Our approach to high energy interactions in QCD is based on the BFKL Pomeron [3] and on 
the reggeon-like diagram technique which takes into account the interactions of BFKL Pomerons 
01 El Ej- It is well known El E] that the Pomeron diagrams technique can be re- written as a 
Markov processes ^0] for the probability of finding a given number of Pomerons at fixed rapidity 
Y. 

However, the colour dipole model provides an alternative approach [H], in which we can 
replace the non-physical probability to find several Pomerons 1 , by the probability to find a given 
number of dipoles. It has been shown that everything we had know about high energy interactions 
can be re- written in the leading N c approximation (N c is the number of colours), in probabilistic 
approach [HI [T21 HH| using dipoles and their interactions. In this approach assuming that we neglect 
the correlations inside the target, we have a non- linear evolution equation (Balitsky-Kovchegov 
equation jlHHS], see Ref- [1211111 for alternative attempts). 

The alternative approaches based on the idea of strong gluonic fields [17], or on the Wilson 
loops approach [Ti] . lead to a different type of theory which is, at first sight, not related to 
Pomeron interactions. For example the JIMWLK equation ^H] appears to be quite different 
to the Balitsky-Kovchegov equation. However, it has been shown that the Balitsky- JIMWLK 
approach, is the same as the method based on the dipole model in the leading N c approximation 
(see Ref. [TTH !2Uj). The further development of the Balitsky - JIMWLK formalism, as well as 
a deeper understanding the interrelations of this formalism with the old reggeon - like technique 
based on the BFKL Pomerons and their interactions, is under investigation (201 El 122] and we 
hope for more progress in this direction. Although, the Balitsky-Kovchegov equation is widely 
used for phenomenology, it has a very restricted region of applicability [221 [2H EH] , since it has 
been derived assuming that only the process of the decay of one dipole into two dipoles is essential 
at high energy. As was shown by Iancu and Mueller [23] (see also Refs. [21 [23 ESI 123 123 EH|), a 
theory neglecting the merging of two dipoles, is not only incorrect, but that the corrections from 
the process of annihilation of two dipoles into one dipole should be taken into account, to satisfy 
the unitarity constraint in the t-channel of the processes. This, in spite of the fact that there is a 
l/Af? suppression in the framework of the dipole approach. 

A number of papers has been circulated recently |2*H1 E7J [2H1 1213 EI] m which the key problem 
of how to include the annihilation of dipoles without losing the probabilistic interpretation of the 
approach is discussed. 

In this paper we find the asymptotic solution to the high energy amplitude using the approach 
of Ref. [2H]- The paper is organized as follows. In the next introductory section we discuss the 
dipole approach in 1/N% approximation, and consider the main equations which we wish to solve 
in this paper. We also discuss the statistical interpretation of these equations which play a major 
role in our method of searching for the solution. The section 3 is devoted to solving the simplified 
toy model where the interaction does not depend on the size of interacting dipoles. We find here 

1 We call this probability non-physical since we cannot suggest an experiment in which we could measure this 
observable. 
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that the large parameters of the problem allow us to use the semi-classical approach for analyzing 
the solution. In section 4 we find the asymptotic behaviour of the high energy amplitude in the 
dipole approach including Pomeron loops. We discuss how the amplitude high energy approaches 
this solution. In the short last section we summarize our results. 



2 Generating functional and statistical approach 
2.1 Dipole model in 1/N* approximation 

The first question that we need to consider is, can we use the dipole approach for calculating I/AT? 
correction. At first sight the answer is negative . Indeed, Mueller and Chen showed in Ref. [SHj 
that the term of order ag/N 2 cannot be rewritten as dipole interaction (see Ref. [20]; i n which 
this result is demonstrated in the framework of the JIMWLK approach). We consider two dipole 
rescattering given by Fig. ^ In the leading order in N c , the gluon can be emitted by a quark (or 
antiquark) of one dipole, and should be absorbed by an antiquark (quark) of the same dipole, as 
is shown in Fig. ^a. In the next to the leading order in N c the interaction of two dipoles has the 
form EDI C2D| : 

-^1 J d 2 zK{x x ,yx,x 2 ,y 2 ,z) (T(x 1 ,y 1 x 2 ,y 2 )) (2.1) 
where T(xi, y\ x 2 , y 2 ) is the interaction operator for two dipoles, as = N c as/n and 

K (xi,yi; x 2 ,y 2 ; z) = — — - — — x — — - — — = 2.2 

V (ari - z) 1 {yi-z) 2 J V ( x 2 - A (y 2 - z) 2 J 

= ]z (K(x 2 ,yi] z) + K(x 1} y 2 ]z) - K{x u x 2l z) - K^y^z)) (2.3) 



where 



K(x,y;z) = ( ! ! )2 (2.4) 
(x — z) 1 (y — z) z 



Eq. ()2.3|) has been discussed in Ref. [SI] in the context of the odderon structure in the 
JIMWLK -approach. The first two terms in Eq. ()2.3|) describe the emission of a gluon by two 
dipoles x 2 — yi and X\ —y 2 (see Fig. [TJb) . The configuration in which quark x\ (antiquark yi) and 
antiquark y 2 (quark x 2 ) creates a colorless pair is certainly suppressed by factor l/N 2 . However, 
after being created dipoles x 2 — y\ and X\ — y 2 , will interact as two dipoles in the leading N c 
approximation (compare Fig. ^-a and Fig. HJ-b). Indeed, the topology of this term is the same as 
two cylinders 2 . Therefore for the first two terms we can replace (T(xi,yxx 2 ,y 2 )) by the product 
{T{x u y 2 )) (r(£ 2 ,&)). 

2 The diagrams describing a dipole target interaction in the leading order in N c has a cylinder topology as was 
shown in Ref. |32j . This means that all diagrams for the BFKL Pomeron exchange can be drawn in the cylindric 
surface. 
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The second two terms in Eq. (J2.3)) stem from the possibility of two quarks rescattering with 
a suppression of This interaction leads to a quite different topology (see Fig. [TJ-c), which 

cannot be treated as two independent parton showers. This new configuration with non-cylindric 
topology should be treated separately using the so called BKP equation jHH] , and it is called 
multi-reggeon Pomeron. Such Pomerons have been studied long ago (see Refs. |3H ESj) and to 
the best of our knowledge, the intercepts of these n-reggeon Pomerons turn out to be smaller than 
the intercept of n cylindrical configurations (in our case n = 2). 

Therefore, we conclude that we can use the dipole model to even calculate l/iV| corrections. 
This analysis is supported by direct calculation in Ref. [B], where it is shown that the triple 
BFKL Pomeron vertex calculated by Bartels in Ref. [3] generates 1/N% corrections which can be 
rewritten as dipole interactions. 



leading N c order 

Xi 



1/N C corrections 



x 2 



C 
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1/N C corrections 

x, _ 




a) 



b) 



c) 



Figure 1: Two dipole rescattering in leading N c order (Fig^a) and 1/N% corrections (Fig.^b and 
Fig. \^-c).Fig. \^b shows the 1/N% corrections due to re-grouping of the dipoles. This contribution 
has the same topology as the leading N c order diagram (see Fig. QJ-a. Fig. Q]-c corresponds to 
interaction of two quarks from the different dipoles. This contribution has a different topology 
than Fig. \J[a, and has to be assigned to so called multi-reggeon Pomerons (see for example Ref. 



2.2 Main equations 



In this section we will discuss the main equations of Ref. [28J paying our attention to their 
statistical interpretations. 

In Ref. ^T] the generating functional which characterizes the system of interacting dipoles was 
introduced, it is is defined as 



Z(Y - Y ; [u]) = (2.5) 

r n 

/ P n (Y -Y ;xi,y 1 \x 2 ,y2] ■■■]Xi,y i ]...;x n ,y n ) Y[ u{xi,Ui) d 2 x { d 2 y { 

— i •> „• 1 



n=l 



i=l 



where u(xi,yi) = Ui is an arbitrary function of Xj and y^. P n is a probability density to find n 
dipoles with the size Xi — y^ and with impact parameter (x« + yi)/2. Directly from the physical 
meaning of P n and definition in Eq. ()2.5|) it follows that the functional (Eq. ()2.5|) ) obeys the 
condition 

Z{Y - Y ; [u = l]) = 1. (2.6) 

The physical meaning of (Eq. ()2.(jj) ) is that the sum over all probabilities is one. 

The functional Z has a very direct analogy in the statistical approach: i.e. the characteristic ( 
generating ) function 0(s) in Ref. [TO] - For P n we have a typical birth-death equation which can 
be written in the form: 

9 P-n ( • • • ; x% j y% ) • • • ) 
&Y 

= -E r i-2 (g) {P n {...,Xi, yi ;... )- P n ^(..., Xl , yi ;...)) (2.7) 

i 

- J2 r 2-i • • > V» ■■■) ~ p n+i(- ■ ■ ; Xi, yr, . . . ; . . .)) 

- E F 2^3 <S> ( p n(- ■■iX i ,y i ;...\ x k , y k ; . . .) - P„_i(. . . , ar 4 , y t ; . . . ; x fe , . . .)) 

In Eq. ()2.7|) r!^ 2 is the vertex for the process of decay of one dipole with size x — y, into two 
dipoles with sizes x — z and y — z. This vertex is well known and it is equal to 

= a s K(x,y;z) (2.8) 

The vertex for the process of transition of two dipoles with sizes X\ — y\ and x 2 — 2/2 into three 
dipoles with sizes x\ — y% % z — y\ and x 2 — z has been discussed in Ref. |2B] and has the form 

r 2 ^s = ^ (K(x 2 ,fcz) + K{x l: y 2 -z)) (2.9) 

c 

We will discuss the 2 — > 1 vertex we will discuss below. 

denotes all necessary integrations (see more detailed form in Ref. [2E])- One can see that 
each line in Eq. ()2.7|) gives a balance of the death of a particular dipole ( the first term in each 
line which has a minus sign), and of the birth of two or three dipoles ( the second term in each 
line which gives a positive contribution . Eq. (J2.7|) is a typical equation for the Markov's process 
( Markov's chain) [TUj . 

The last term in Eq. (|2.7|) has the following form: 

T 2 ^3 (g) {P n {- ■ ■ , Xi, y^, ■ ■ ■ ; x k , y k ; . . .) - P n _i(. . . , x h y { ; . . . ; x k , y k ; . . .)) = (2.10) 

= J d 2 Xid 2 x k d 2 yid 2 y k d 2 z (K(xi,y k ;z) + K(x k ,yi,z)) 
\Pn\- • • 1 Xi, yi, ■ ■ ■ , x k , y k , . . .) i 3 ra _i(. . . , Xi, yi, . . . , x k , y k , . . .)) 
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Multiplying Eq. (|2.7|) by the product n£=i u i an d integrating over all Xj and yj, we obtain the 
following linear equation for the generating functional: 

dZju] 
dY 



with 



= x[u]Z[u) (2.11) 

XM = (2.12) 

- J d 4 qd\d% \T x ^ 2 {q^ qx + q 2 ) (-u(q) + u( gi ) u(q 2 ) ) - (2.13) 

-r M + ft - ff )(«(gO «(*)-«(«)) L( 9l )L(J - (2-14) 



2 3 



~ / II ^ f[ rf 2 ^ r 2 ^ 3 ((a*, + (x 2 , y 2 ) -> (xi,y 2 ) + (z, yi) + (x 2 , z)) 



t=l i=l 



1 <5 2 

(«(a?i,Jfe) - 1) (u{x2,yi) - u(z,yi)u{x 2 ,z)) ~ x , . x( : (2.15) 

2 du(gi) dw(g 2 ) 

where we denote d A qi = d 2 Xid 2 yi. The functional derivative with respect to u(q) = u(x,y), plays 
the role of an annihilation operator for a dipole of the size r = x — y, at the impact parameter 
b = |(x + y). The multiplication by u(x,y) corresponds to a creation operator for this dipole. 
Recall that d 4 qi stands for d 2 X{ d 2 yi. 

In Eq. ()2.12|) we subtracted the term that corresponds to the 2 — > 3 transition at u(xi, y 2 ) = 1. 
Indeed, at u(xi,?/ 2 ) = 1 this transition describes the decay of the colour dipole (x 2 ,yi) into two 
dipoles at any size of the dipole {xi,y 2 ) over which we integrate. Such a process has been taken 
into account in the first term of Eq. (|2.12j) which accounts for 1 — > 2 decays of all possible dipoles. 

3 

Eq. (|2.11|) is a typical diffusion equation or Fokker-Planck equation [TU], with the diffusion 
coefficient which depends on u. This is the master equation of our approach, and the goal of 
this paper is to find the asymptotic solution to this equation. In spite of the fact that this is a 
functional equation we intuitively feel; that this equation could be useful since we can develop 
a direct method for its solution and, on the other hand, there exist many studies of such an 
equation in the literature ( see for example Ref. JUj) as well as some physical realizations in 
statistical physics. The intimate relation between the Fokker-Planck equation, and the high energy 
asymptotic was first pointed out by Weigert [SH] in JIMWLK approach, and has been discussed 
in Refs. jSH 123 HZ1 • 

It should be stressed that in the case of leading N c approximation, the master equation has 
only the first term with one functional derivative and, therefore, the Fokker-Planck equation 
degenerates to a Liouville's equation and describes the deterministic process, rather than stochastic 

3 In the first version of this paper as well as in the first version of Ref. [2B1 we made a mistake of forgetting 
about this term. In doing so, we incorrectly generated the two Pomerons to one Pomeron transition. We thank all 
our colleagues whose criticism helped us to find a correct form for 2 — ► 3 transition in Eq. (|2.11|l . 
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one which the Fokker-Planck equation does. The solution to the Liouville's equation is completely 
defined by the initial condition at Y = 0, and all correlations between dipoles are determined by 
the correlations at Y = 0. As has been shown ^21 EH ED that only by assuming that there are 
no correlations between dipoles at Y = 0, we can replace the general Liouville equation by its 
simplified version, namely, by non-linear Balitsky-Kovchegov equation f3J E] • 



2.3 r(2->l) 

For further use we need more detailed information on the vertex, for merging of two dipoles in 
one dipole. 

As was shown in Refs. [2*71 128) 12*9"] the vertex r (2 — > 1) can be found from the integral equation 
which has the form 

d 2 xd 2 yT 2 _ rl (x 1 ,y 1 + x 2 ,y 2 -»• x,y) >y BA (x,y;x',y') = (2.16) 

d 2 x[, d 2 y[ , d 2 x' 2 , d 2 y' 2 T^ 2 (a/, y' -> x[,y[ + x' 2 , y 2 ) >y BA (x 1: y x ; x[, y[) >y BA (x 2 , y 2 ; x' 2 , y' 2 ) 

where •y BA is a dipole-dipole elastic scattering amplitude in the Born approximation, which is 
equal [31 EE] 

BA, i t\ «! _JL_ i 2 ( (x - x') 2 (y - y 1 )^ 
7 &V> x >y) = 2iVf 16^ ln \ (x - y'Y (y - x>y ) ^ 

Eq. 1)2.16]) is the basic equation from which the vertex T 2 ^x can be extracted. To achieve 
this we need to invert Eq. ()2.16j) . by acting on both sides of it by an operator which is inverse 
to '-/ba in operator sense. Fortunately, this operator is known to be a product of two Laplacians 

ma 1231211 mi: 

2N 2 r 

r 2 ^x{x 1 ,y 1 + x 2 ,y 2 -*x,y) = - z fA x A y d 2 x[,d 2 y[ ,d 2 x' 2 ,d 2 y' 2 x (2.18) 

a s J 

r^ 2 (x',y' -> x' v y[ + x' 2 ,y 2 ) >y BA (x u y x \ x[, y[) >y BA (x 2 , y 2 ; x 2 , y' 2 ) 

The exact evaluation of Eq. 1)2.18)1 is done in Ref . [2H] , but for further presentation in this paper 
we will need the vertex T(2 — > 1) only in the form of Eq. 1)2.18)1 . 

2.4 Scattering amplitude and its statistical interpretation 

As was shown in Refs. [i"""] I13j) the scattering amplitude is defined as a functional 

oo f^n n g 



\Ui=l 



N (Y; [ 7 ,]) = / £ IT **i *Vi 7fo, Vt) ^ (Y, M) 

J n=1 n. i=1 ou % 

„ oo n 

= ~ J El" 1 )™ IT d 2 x l d 2 ya(x l ,y l )p n (r 1 ,bu...r n ,b n ;Y ) (2.19) 



n=l i=l 
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where 7 denotes an arbitrary function which should be specified from the initial condition at 
Y — 0. To calculate the amplitude we need to replace each term ri?=i l{ x h Hi) by function 
7„(xi, yx] . . . ; x n , y n ) which characterizes the amplitude at low energies for simultaneous scattering 
of n dipoles off the target. 

From Eq. ()2.19j) one can see that 

1 - 5 

p n {xi,yi;...;x n ,y n ) = -t J]- — Z (Y; [ui]) \ Ui=1 (2.20) 

n\ f-J[ d Ui 

As was shown by Iancu and Mueller t-channel unitarity plays an important if not crucial 
role in low x physics (see also Refs. (2HI2S1)- In the context of this paper it should be noted that 
t-channel unitarity as a non-linear relation for the amplitude, is able to determine the unknown 
parameters in the asymptotic solution, t -channel unitarity for dipole-dipole scattering can be 
written in the form (see Fig. EJ) 

N(x,y;x',y';Y) = (2.21) 

00 „ n 

X) / II d2%i d2 y* ^ p ( x ' v> ■ ■ • ; x *> • • ■ v™ Y ~ Y ') v'\ ■ ■ ■ ; x ^ !/<;••■ x n, y n ; y') 

n=l J i=l 

where N p and N l denote the amplitude of the projectile and target, respectively. Actually, they 
are not exactly the amplitude. Indeed, accordingly Eq. (j2.21j) their dimension should be 1/x 2 while 
the amplitude iV is dimensionless. On the other hand, we know that the unitarity constraint has 
a form Im N = V] n N{2 — > n) n*(2 — > n) illustrated in Fig. |21 However, the amplitude in this 
relation should be taken in the momentum representation, while we here consider the amplitude 
in the coordinate representation. The difference is clear from Eq. (j2.19|) where we have that the 
amplitude of interaction at low energy 7(2^, yi) has been taken into account in the definition of N. 
On the other hand in the unitarity constraints, these amplitudes should only be included once for 
both amplitudes (see Fig. EJ). The general way to do this is to re- write the unitarity constraints 
through the function p n using Eq. ()2.2()j) . 

It leads us to the following form of the t-channel unitarity constraint: 

N (x,y,;x r ,y';Y) = 



= £ (-l)»n! J II d 2 x> S Vl J] d 2 x[ d 2 y> 7^^; 4 2/0 (2-22) 

n=l i=l i=l 

P^n ( x i V\ x i-> 2/i 5 • • • x ii Ui] ■ ■ ■ j x n yni Y — Y ) p n (x ,y ] x 1i y^, . . . x i , y^ . . . ; x n , y n ; Y ) 

The factor n\ appears in Eq. ()2.22j) is due to the fact that each dipole with rapidity Y' from the 
target can interact with any dipole from the projectile (see Fig- EJ) • 

The unitarity constraint itself shows that the high energy amplitude could be described as 
a Markov process. Indeed, this constraint claims that the amplitude at later time ( Y + AY) 
is determined entirely by the knowledge of the amplitude at the most recent time (Y) since 
N(Y + Ay) = N(AY) N(Y) from t-channel unitarity. 



N(Y) = £ ] 
n 1 




N(Y-Y') 



N*(Y') 



Figure 2: The unitarity constraint 
in t-channel. 



The relation between Z and N is very simple |15j . namely 

N(Y,[ 7i }) = 1 - Z(Y,[1- 7j ]) 



(2.23) 



Therefore, the amplitude N is closely related to the characteristic function in the statistical 
approach. However, it is well known that it is better to use the cumulant generating function 
which is the logarithm of the characteristic function. In our case, we introduce the cumulant 
generating functional, namely, 



oo „ n 

$(X;N) = foZ(Y;[ji\) = / II d 2 Xid 2 yai i^hVi) p n (• • • • • •) 

~, 1 J A 1 



(2.24) 



n=l i=l 



where Ui — 1 — 7*. Eq. (|2.24|) is the definition for functions p n which are related to cumulants. For 
example p 2 is related to (|T(r 1 )T[r 2 )|) — (|T(r 1 ))(|T(r 2 )) rather than to (|T(rx)T[r 2 )|) as P n does. 
T is the scattering operator. 

The advantage of using $ is the following: (i) if there are no correlations between dipoles it 
is necessary and sufficient to keep only the first term in the series of Eq. ()2.24|) ; (ii) to take into 
account the two dipole correlations we need to keep the two first terms in Eq. (J2.24)) ; and (iii) the 
n-term in the series of Eq. (J2.24|) describes the correlations between n-dipoles. 

Substituting Eq. (|2.24|) in Eq. (|2.11|) we obtain the equation for $ 



d$(Y; N) 
dY 



J d 4 qd%d% T(q ^q 1 + q 2 ) ( 7 (<?) - 7 (?i) lifb) ) $ (Y; N) + (2.25) 



+ \j AAl^lXgi + q 2 -> q) (7(g) - 7(^1)7(^2)) X 



6 2 



5 7 (g 2 



Sj(qi) 5<y(q 2 
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Figure 3: An example of 
enhanced diagrams that 
contribute to the unitarity 
constraint in t-channel. 
Zigzag lines denote the 
BFKL Pomerons. This 
particular set of diagrams 
can be summed by the 
unitarity in t channel 
and corresponds to Iancu- 
Mueller approach 



+ 9 / Y[d 4 q i d 2 zT( y q 1 + q 2 -> (x u y 2 ) + (x 2 , z) + (z, y^) j(x x , y 2 ) (l(x 2 ,y 1 ) - j(x 2 , z) 7(3, y%) ) x 

(^* (y;W) ^* (r;W) + ras)* (r;W) ) 

where u(q) = 1 — 7(g). 



3 Asymptotic solution in the toy-model 

3.1 General description 

We start to solve the master equation (see Eq. (j2.11|) and Eq. ()2.25|) ) by considering a simple 
toy model in which we assume that interaction does not depend on the size of dipoles (see Refs. 
(TTJElEl! f° r details). For this model the master functional equation ( see Eq. (j2.11j) ) degenerates 
into an ordinary equation in partial derivatives, namely 

dZ 

w = ( 3 - 26 > 

= -r(l -»• 2) u(l - u) ^- + r(2->l) u ( l-u) ^-=- + r(2 -> 3)u (1 - u) 2 ^- 



du (du) 2 " (du) 
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To obtain the scattering amplitude we need to replace 7 in Eq. ([2.23)1 . by the amplitude of 
interaction of the dipole with the target. For N Eq. (J3.26)) can be rewritten in the form: 



dN 

w 



(3.27) 



dN 



r(i - 2) 7(1 - 7 ) + r(2 ^ 1) 7(1 - 7 ) 



o 2 n 



if 7 is small we can reduce Eq. f!3.27|) to the simple equation 



+ r(2^3) 7 2 (l-7) 



o 2 n 



dN 
~dY 



rn 



2)7 



dN 



(3.28) 



(3.29) 



which has the solution 

N = 7 e r ^ 2 ) y 

which describes the exchange of the Pomeron with the intercept T(l — > 2). 

All other terms in Eq. 1)3.27)1 have a very simple physical meaning, describing Pomeron inter- 
actions (see Fig. HJ). 



P = e 



ra->2)Y 




= - r(2->3) 



p 
p 



/ 



= -r(i->2) 



p p 

X= r(2->3)-r(2->i) 
sP 



Y 



=r(2->i) 



Figure 4: Pomeron interactions described by Eq. \3.21 ). 



The probabilistic approach that we develop in this paper has an advantage that we can apply 
the well developed formalism of the statistical physics 1(> to the partonic system (see Refs. 
(HI E3 12H1 121] where statistical approach has been applied to our problems). It is reasonable to do 
this if we are dealing with a system with a large number of dipoles at high energy. We will try to 
illustrate this point below, but first let us describe the strategy that we will follow in our search 
for a solution. 



3.2 The strategy of our search for the solution 

Generally speaking the solution to Eq. (J3.2(-)j) or to Eq. ()3.27j) can be written in the form 

Z{Y;u) = Z(oo;u) + AZ(Y,u) (3.30) 
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where asymptotic solutions Z(oo;u) or N(oo;u), are solutions of Eq. (|3.26() or Eq. (|3.27() with 
the l.h.s. equal to zero, namely, it means that asymptotic form of Z can be determined from the 
equation 

- r(l->2)«(l-«)2^) + (3.31) 

+ r{2 ^ 1)u(1 _ u) ll^p. + r ( 2^3) M (l- M r^^)=0 

If we prove that the solution AZ(Y,u) which satisfies the full Eq. (13.26(1 decreases at high energy 
( at Y — > oo) the asymptotic solution will give us the behaviour of the scattering ampllitude at 
high energy. 

Eq. (|3.26() and Eq. (|3.27() are diffusion equations, in which the diffusion coefficient is a function 
of u or 7. Let us solve the simplified problem replacing Eq. (|3.26() by the equation with constant 
coefficient, namely, 

dZ _dZ ld 2 Z 
dy du k du 2 

where y = T(l — > 2) Y and k is defined in Eq. (I3.3fij) . 
Then 

Z(oo,«) = e K(u - 1} (3.33) 

in which we fixed all constants, assuming that all other solutions give small contributions at large 
values of Y and Z(oo; u = 1) = 1. 

The general solution for Eq. ()3.32() is known, namely, it is equal to 

AZ(Y;u) = j-V^expf-/ ^ (3.34) 

Eq. (|3.34() is the Green function of Eq. ()3.32() . and using it one can construct the solution for 
any initial condition. This solution vanishes at large y. This fact shows that Eq. ()3.33() is the 
asymptotic solution at high energy. 

Therefore, the strategy of searching for the solution to our problem consists of two steps: 

1. Finding the asymptotic solution as a solution to Eq. ()2.11j) with zero l.h.s.; 

2. Searching for a solution to the full equation, but assuming that this solution will be small 
at high energies (Y — ► 00). 

We also have to show that the solution, which we find, satisfies the initial condition which we take 
Z{Y = 0,u) = u. 

In the following presentation we consider separately two cases T(2 — ► 3) = and T(2 — >• 1) = 
which lead us to further understanding of the asymptotic behaviour of the scattering amplitude. 
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3.3 Solution for r(2 — ► 3) = 
3.3.1 Asymptotic solution 

Considering a simple case with T(2 — > 3) = we obtain a solution for Z(oo; u) from Eq. (j3.31|) 

Z(u; Y -»• oo) = 1 - 5 + Be^ 1 ) (3.35) 

where 



T(l -> 2) 2 AT 2 HI — 

It should be stressed that in the toy model we keep the same order of magnitude for the vertices 
as in the full QCD approach. This explains all numerical factors in Eq. ()3.36|) . To avoid any 
confusion we would like to draw the reader attention to the fact that we will use k and k in 
further discussions to denote the combination of as and N c but not the ratios of the vertices 
which are functions of coordinates in the general case. These combinations give the order of 
magnitude for the ratios of vertices as far as the a$ and N c factors are concerned. 

Eq. (J3.35|) satisfies the initial condition that Z(u = 1;Y) = 1, but the coefficient B remains 
undetermined, since we cannot use the initial condition at Y — Yq. We will show below that the 
unitarity constraint in t-channel will determine B — 1. Expanding Eq. ()3.35|) we obtain: 

lim P (Y) = e" K ; lim P n (Y) = — e~ K ; (3.37) 

Y^oo Y— >oo n\ 

We obtained the Poisson distribution with the average number of dipoles (n) = k ^> 1. Hence 
we will attempt to discuss the evolution of our cascade using statistical methods 4 . 

To find the value of B in Eq. (|3.37|) we can use the unitarity constraint in t-channel, which 
can be obtained from Eq. (|2.22jl . and which for the simple model reads as EHl I2U IT3] 

OO 1 

N(Y) = Y,{-l) n ^—p^{Y-y)fi{y) (3.38) 

n=l K 

where p both for the projectile and the target is defined by Eq. ()2.21|) . Eq. ()3.38j) corresponds to 
a simplification that in the toy model r y BA = ^ 8(xi — x^) 5(yt — y[). 

Rewriting Eq. ()3.38|) for Y — y ^> 1 and y ^> 1 and taking into account Eq. ()2.23|) we have 

OO 1 

iV(oo) = £ (_)»„! ^(oo) p^(oo) (3.39) 

n=l K 

Therefore, we see that B = 1 since for large k Eq. (|3.39|) reads as 



B = B 



2 



'in this subsection we follow Ref. 9 in which a wide class of models for interacting Pomerons is considered. 
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Substituting p^(oo) = «"/n! and p^(oo) = we obtain 

iV(oo) = 1 - e"^ = i _ e -« (3.40) 

where k p and are related to the projectile and target, and they are equal k v = k 1 = k in 
our case. The unitarity constraint determines the value of u in the solution of Eq. ()3.35|) which 
is relevant to the asymptotic behaviour of the amplitude. We can translate the result of using 
unitarity given by Eq. (|3.4(Jj) . as the u = contribution to the asymptotic behaviour of the 
amplitude, if we calculate it from Eq. (|3.35|) . 



3.3.2 Search for energy dependent solution 

Eq. (|3.35|) leads us to use the Poisson representation 0: 

P{Y) = J daF{a,Y)(^e~ a ^ = (P n (a)) (3.41) 

where < . . . > stands for averaging with weight F(a, Y). 
The generating function Z(Y; u) can be written as 

Z(Y-u) = J da F(a,Y)e {u - 1)a (3.42) 

Our master equation (see Eq. (|3.26j) ) can be rewritten as 

d F(a Y) d Id 2 

—^y- 2 = -Q- a (A{a)F(a,Y)) + (B(a)F(a,Y)) (3.43) 

where 

A(a) = r(l -»• 2) a - T(2 -> 1) a 2 ; (3.44) 
B(a) = 2A(a); (3.45) 

Eq. f!3.43|) is an equation of Fokker-Planck type ( see Refs. EH1 12H1 12Z| with positive diffusion 
coefficient B(a) at least for a < k. The initial condition of Eq. (J2.6|) gives the normalization of 
F(a, Y), namely, 

J daF(a,Y) = 1 (3.46) 

For positive B(a) the distribution function F(a,Y) is positive with the normalization of 
Eq. ()3.4(ij) . and therefore can be considered as a probability distribution [TQJ . Hence, the en- 
ergy dependence of the scattering amplitude can be given by averaging of Eq. ()3.40|) . namely, 

N(Y) = (3.47) 
= ((iV (a ' a )) = J daF(a,Y-y) J daF(a,y) (l - e~ a ") 
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Eq. ([3.43)1 is equivalent to the differential stochastic equation [TU] 

da = A(a) + ^B(a)dW(Y) (3.48) 

where dW(Y) is a stochastic differential for the Wiener process. Using Eq. (J3.48)) . together with 
the large number of dipoles involved in the process is the reason we hope to solve the whole 
problem using the statistical approaches (see Refs. 123 Hi] ) • However, it is too early to judge 
how successful this approach will be. 

As one can see from Eq. ()3.48|) a grows until it reaches the zero of B(a) and it becomes 
frozen at this value at large Y . Therefore, we can solve Eq. (|3.43|) at large energy assuming that 
k - a < /t. In this limit A (a) can be replaced by A (a) = k — a = £ and Eq. (|3.43j) reduces to 
the form 

9F(a) , _ . dFioc) / _d 2 F(a) . „ . 

where y = T(l —>■ 2) Y with the solution 

F{^y) = I ^e« y F(a,u); (3.50) 

and 

F{ot,u) = 0(^)1^(1 + ^,2,0 (3.51) 

where 1F1 is the confluent hypergeometric function of the first kind. We should choose the function 
4>(lj) from the normalization condition of Eq. ()3.46|) and, finally, 

F{a,L>) = ^F, (1 + ^,2,0 (3.52) 

To find the asymptotic behaviour of this solution at Y ^> 1, reconsider Eq. (|3.5U|) which has the 
form 

F(a, Y) = fipL^ ra-a) v ±Fi (1 + 2 (3 . 53) 

We can use the integral representation for iF% (see 9.211(1) in Ref. [SH]) and return to the Poisson 
representation of Eq. ([3.42)1 . Doing so, the solution has the following form 

Z(u, Y) = e K(n - 1} * (y + In (j^)) (3-54) 

where t = 2(u — 1) and the arbitrary function \1/ should be determined from the initial condition. 
One can see that from Z(u = 1;Y) = 1 this function is equal to 1. 

This solution corresponds to the asymptotic solution which does not depend on Y (see Eq. (|3.35|1 ). 
We have to consider a different region of a to find out how our system approaches the asymptotic 
regime given by Eq. (|3.35|) . Assuming that a k we can model A and B in Eq. (|3.43J) by 
A(a) = a. In this case Eq. (13.43)1 has a simple form 

d F(a Y) d d 2 

= (aF{aX)) + _ (aF{aX)) (3.55) 
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Returning to the variable t = u — 1 (see Eq. ()3.42|) ) we obtain the equation for Z 

d -§£l = -2Z( t ,y)-t<H-*>^P (3-56) 

which has the solution 

Z(.,y)= <i + £ S ( y + lB (<i±4)) (3.57) 

The solution of the master equation (see Eq. ()3.26j) ) can be written as 

Z(t,Y) = e«> + (i±^ S ( y + ,„((i±^)) (3.58) 

Function 5 is arbitrary function and we determine it from the initial condition that Z(u,Y = 
0) = u. Finally, the solution has the form 

Z (t ,Y) = e- + e-<-^- ^ ——LtL—— _ exp{K _j£^_ } ) (3.59) 

Therefore, the key difference between the case with only emission of dipoles (only T(l — > 
2) 7^ 0), and the case when we take into account the annihilation of dipoles T(l — > 2) ^ 
and r(2 — > 1) 7^ 0, is the fact that there exists an asymptotic solution which depends on u (see 
Eq. (|3.35Jl ). Eq. (|3.59|) shows that the energy dependent solution in the wide range of t or u satisfy 
the initial condition Z(Y = 0; it) = u and it decreases at Y — ► 00. Therefore, Eq. (|3.4(Jj) gives the 
asymptotic solution to our problem. 



3.4 Solution for T(2 — ► 1) = 
3.4.1 Asymptotic solution 

For the case r(2 — > 1) = the system shows quite different behaviour, namely, 

Z(u; Y -> 00) = 1 + — B — ((1 - uf- k - (1 - no) 1 ^) (3.60) 

1 — K v ' 

where k is given by Eq. ()3.36|) . Formally speaking we cannot satisfy the boundary condition 
Z(u = 1;Y — > 00). The reason for this is obvious, since we cannot take u close to 1 and neglect 
r (2 — > 1) term in Eq. f)3.26j) . We can do this only for u > u Q ~ 1 — k/k oc a|. 

One can see that Eq. ()3.60|) gives 

< n >= dZ(u; Y 00)/ du\ u=l = B(l - u )-~ K B {R/n)' k > 1 
Therefore, once again we hope that the statistical approach can work in this system. 
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Eq. (J3.6U|) is the asymptotic solution in this case. Once more we should use the unitarity 
constraint to fixed the parameter B. Repeating all calculations that led to Eq. (|3.40j) we obtain 

~ 1 T 2 (k + n) ( 1 \ n , 

= 2 Fo\k,k;; — — J < 



< 



where T(x) is the Euler gamma function, 2 F is the generalized hypergeo metric function 13*!?] and 
Ei(x) = j^oo \dt is exponential integral. Considering 1 — u ~ — k/k we see that 

iV(oo) < 1 - exp + exp (^) Ei (-^) (3.62) 



3.4.2 Solution for AZ(Y; u) 

As has been shown in Eq. (|3.6U|) . the asymptotic distribution at large values of Y is not the 
Poissonian type. If we try to introduce the distribution function F(a, Y), the equation for it has 
the form 

dF{a,Y) 



dY 



(3.63) 



f) 1 r) 2 <9 3 

-- (A(a)F(a,Y)) + -_ (B(a)F(a,Y)) - — (C(a)F(a,Y 

where Y = T(l -> 2)Y and 

Ala) = a - - a 2 ; (3.64) 

1 1 \ 2 



B(a) = 2 a -2 a ; (3.65) 

\ K K J 

C{a) = la 2 ; (3.66) 

K 

B(ol) does not have a zero and, therefore, we expect that the asymptotic behaviour will be related 
to the large values of a. For the master equation (see Eq. ()3.26j) ) it means that u — > 1 will be 
essential. In the simple case T(2 — > 1) = and u — ► 1(1 — u <C 1 the master equation degenerates 
to 

dZ(u,Y) 



T(1^2)dY 

u(1 u) dZ(u,Y) + (1-n) 2 g 2 Z(n,7) ^ l ^CC.r) + 



(3.67) 



du k du 2 k dC, k d(, 



2 
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where £ = ln(l — u). The solution of this equation is simple, namely, 

ra+ioo d\ 



Z(C,Y) = P lC ° ^±Z(X)e(^ x + ^ x2 ) r{1 ^ Y + x < (3.68) 

J a— ice 2 7T 2 



where Z(X) = 1/(A — 1) using the initial condition that Z(u,Y = 0) = u. 
The solution which satisfies our initial condition is equal to 



Z((,Y) = erf I I (3.69) 



2^r(i 2)y/« 

where rj = (1 + 4) T(l — > 2) Y + £ and erf(x) is the error function given by erf(x) = 
(2/0F) JJ exp(-/) dt. 

At K = Eq. ()3.69|) leads to Z((, Y = 0) = 1 ps it, since we assume that w is close to unity. 

At K — > oo we have 



Zfey oo) - 1 - . illl^lZ exp f- (3.70) 
W ' ; y nkri P V 4F ( 1 2 ) */ 

Therefore, we have obtained a solution and the only problem, that remains, is the behaviour 
of the asymptotic solution at u — > 1. We need to investigate the region of u — > 1, but we will 
do this using a new approximation, in which we use the large parameters of our approach ( see 

Eq. (nrnj)). 

3.5 Semi-classical approach for large k and k 

Considering the solution of Eq. (|3.59|) . we notice that the function $ (see Eq. ()2.24|) and 
Eq. (|2.25jl ) for this solution is large and it is proportional to k. This observation triggers a search 
for a semi-classical solution assuming that <3>^ >> $ 7i7 where $ 7 = dPhi/d^. 

For r(2 — > 3) = in the toy model Eq. (|2.25|) has the form 

+ (,n, 

Fist, we see that the asymptotic solution $(0057) = —kj is the same as that we obtained in 
Eq. ()3.35|) and Eq. ()3.59|) . We try to find the solution to Eq. (j3.71J) for all values of Y in the form 
®(y;i) = $(00; 7) + A$(3 ; ; 7) considering A$ < $(00; 7). For such a solution Eq. (j3.71j) 
reduces to a linear equation 

^2)._ 7(1 _ 7) ^i2l2) (3.72, 
ay 0*7 

The solution to this equation has the form 

A$(y, 7 ) = H[y+ lnfl— ^]] (3.73) 



7 
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where H is an arbitrary function which should be found from the initial conditions Z(Y = 0, [itj]) = 
u. 

Assuming that A$ -C $(0x0,7) at Y = 0, we can reconstruct the solution given by Eq. ()3.73|) . 
namely, 

Z(~f,Y) = exp($(oo, 7 ) + A$(F, 7 )) = (3.74) 

1-7 , 7e- r(1 - 2)y , 
exp{-/€7 + k- ■ zwt^vF J 



l-7 + 7e- r (^ 2 ) y ' l- 7 + 7 e- r ( 1 ^ 2 ) y - 

This solution differs from the solution given by Eq. (|3.59|h but leads to the same behaviour at 
large values of Y . The difference is obvious since A$ is not small at small values of Y. 

We attempt to find a solution using the same methods as for the case r(2 — ► 1) = 0. 

The asymptotic solution for $ is $(00; u) = —fclnfr). This can be seen directly from 
Eq. (|3.fi()j) . and it also appears as the solution to the following equation 

0. 7(1 _ 7){ « + 1 7 («) 2 } ,3,5, 

For A$(Y» ($(y,u) = $(00; u) + A$(F, u)) we have 

S^m^ 7(1 - 7 )^i2^ (3.76) 
aj^ c/7 

A general solution to this equation is 

A$(F, 7 ) = S (y - In (-^-) j (3.77) 

where S is an arbitrary function which should be determined from the initial conditions Z(Y = 
0; u) = u which translates into initial conditions for A$ as 

A$(F = 0,m) = «ln(7) +ln(l- 7 ) (3.78) 
Finally, we obtain the solution in the form: 

Z(y ;i ) = 1 + - (l - 7 - ie~ y Y~ K (3.79) 



1 — 7 + je 



This is the solution which tends to the asymptotic solution as Y — > 00, at least at small values 
of 7 < 1/R. Note, we have not achieved the correct normalization Z(u = 1,Y) = 1. However, we 
believe this is connected to the problem of ill defined limit of T(2 — > 1), for small 7 < k/k. To 
get a better understanding the situation better we consider the general case in the semi-classical 
approach. 

3.6 Semiclassical approach to a general case 

It is easy to obtain the asymptotic solution 

$(00,7) = -Kln(l + ^7^) (3.80) 
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A<&(y» ($(Y» = $(oo;w) + A$(Y,w)) has the form 

ay c?7 

and the solution to this equation has the same form as Eq. ([3.77)1 . but with a different initial 
condition: 

A$(F = 0;7) = Rhi(l + | 7) + ln(l- 7 ) (3.82) 

Using the arbitrary function S ( see Eq. f[3.77|) ) and from Eq. ()3.82|) we, finally, obtain the 
answer for the generating functional Z 

This solution satisfies all requiments: (i) Z (Y, u = 1) = Z (Y, 7 = 0) = 1; (ii) Z (Y = 0, it) = 
it = 1 — 7j and at Y — > 00 Z ( Y, -u) — > Z(oo, -u). 

The asymptotic solution of Eq. ()3.8()j) leads to the following Z 



^(00,7) = ( 1 + 5 ) (3.84) 



Using Eq. (|2~2T1|) we obtain 

- - 0" 

This p n is the same as for the case of T(2 — > 1) = (see Eq. (|3.6L)jl ) if 1 — Uq is chosen to be equal 
to 1 — u = k/k. It should be stressed that Eq. (|3.83|) and Eq. (|3.85|) lead to the scattering 
amplitude given by Eq. 1)3.61)) . It shows that we correctly guessed the value of 1 — uq. Only for 
1 — it < 1 — u Q do we have to take into account the process of merging of two Pomerons into one 
Pomeron. In other words, we can neglect the 2 — > 1 process if we put the initial condition for the 
generation functional Z(Y,u) at u — uq, namely Z(Y,u — 1 — k/k) — 1. This property resembles 
and even supports the idea of Ref. |45J, that merging processes suppress correlations when these 
correlations are small. 



3.7 Lessons for searching a general solution 

We view this toy model as a training ground to help us find a general solution, and as an aid 
suggesting directions for such a search. We learned several lessons that will be useful: 

1. Our approach in searching for the solution consists of the following steps: 

• Finding the asymptotic solution as the solution to the master equation ( see Eq. ()2.11)) . 
Eq. (I2~231) . Eq. ijsA and Eq. iCPBt ) with zero l.h.s.; 
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• Using the large parameters of our theory given by Eq. (|3.36|) , we can develop the semi- 
classical approach for searching both for the asymptotic solution and for the correction 
to this solution, that provide the form of the generating functional approaching its 
asymptotic value; 

• The corrections to the asymptotic solution decrease at large values of Y, and can be 
found from the Liouville-type linear equation; 

• The important region of u (7) are u — > 1 and 7 — > 0, which should be specified by 
using the unitarity constraint; 

2. The inclusion of T(2 — > 3) is very crucial for the form of solution; 

3. The asymptotic form of the solution differs from the solution to the Balitsky-Kovchegov 
equation: amplitude iV does not approach 1 (N — > 1) at Y — ► 00, but rather N(oo, 7) is a 
function of 7. In particular, it means that we do not expect geometrical scaling 0Hj for the 
solution to the full set of equations. We also do not expect that the scattering amplitude to 
show a black disk behaviour on reaching the value of unity at high energies. We predict the 
gray disc behaviour for this scattering amplitude; 

4. It is easy to find the solution using the equation for the cumulant generation functional 
$ (see Eq. fl2~2lj) and Eq. flZSfl ). 

5. The role of the 2 — > 1 process is very interesting: this process suppresses the small values 
of 7 < k/k in a such way that they can be neglected, supporting the idea of Ref. 

It should be stressed that we found such a solution to the simplified evolution equation (see 
Eq. (J32BD), which has the form of Z(oo,u) + AZ(Y,u) with AZ(Y,u) -> at large Y. It is 
shown that AZ(Y, u) satisfies the initial condition for the generating functional. Therefore, for the 
toy model we prove that our solution is the solution to the evolution equation if we believe that 
we have the only one solution. For the solution to the general evolution equation (see Eq. (j2.11|) 
) we will not be able to show that our solution satisfies the initial condition for the generating 
functional. Therefore, it is possible that could exist a different solution which we overlooked in 
our approach but which will satisfy the initial condition while ours does not. This is the reason, 
why the solution to the evolution equation in the toy model, in which we can check that such 
situation does not occur, is so important for our approach. 

It should be stressed that all new feathures of the solution appear only in the next to leading 
order in 1/N C . In the leading order the JIMWLK-B approach leads to a solution with behaves as 
a black disc. 

At first sight, the gray disc behaviour of the scattering amplitude looks unrealistic. However, 
the situation is just opposite, the gray disc behaviour is rather natural in the parton model (see 
Ref. 50] ) and we have spent a decade to understand how it is possible that we have a black disc 
behaviour for the scattering amplituide in QCD [3T]. The black disc behaviour stems from the 
simple fact that in JIMWLK-B approach the number of 'wee" dipoles of each size increases as 
power of energy due to BFKL emission. It means that the dipoles started to interact with each 
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other even if this interaction is small (proportional to as)- At very high energy the probability to 
find any dipole is equal to unity The Pomeron loops lead to diminishing of the BFKL Pomeron 
intercept and, therefore, to a suppression of the number of the 'wee' partons. Finally, the increase 
of the number of 'wee' dipoles stops before the probability to find a dipole reaches 1. The entire 
picture seems to be close to so called critical Pomeron scenario (see Ref. |S2]): the only theoretical 
model of Pomeron that has been solved. 

Armed with this knowledge we will now to solve the general equations ( see Eq. (j2.11J) and /or 
Eq. (|2~23jl ). 

4 Asymptotic solution (general consideration) 

4.1 Solution for r(2 — > 3)=0 
4.1.1 Solution at Y — > oo 

Our strategy in searching for the solution will be based on the lessons we learnt from the toy 
model. First we try to solve the master equation (see Eq. ()2.25j) ) assuming that 

< —. r—. r (4.86) 



yi)5-f(x 2 , y 2 ) v 8-f(x x ,y x ) 8-f(x 2 ,y 2 ) 

$(oo;[7i]) = -k J d 2 xd 2 y A x A y <j)(x,y; 00) j{x,y) (4.87) 

Substituting Eq. (|4.87|) into Eq. ()2.25|) and putting the l.h.s. of Eq. ()2.25|) to zero, we obtain the 
following equation for asymptotic solution <f>(x, y;Y = 00): 

= J d 2 z K(x, y; z) A x A y <f>(x, y; 00) - A x A y ^J d 2 z K(x, y; z) J d 2 x x d 2 y x (4.88) 
^ BA (x, z\ xi yi) A X1 A yi (f)(xi, yr, 00) J d 2 x 2 d 2 y 2 ^ BA (x, z; x 2 y 2 ) A X2 A w 0(a;i, y x \ 00) 



where K is given by Eq. (|2.4j) and we used Eq. (|2.18|) for the vertex T(2 — > 1). Integrating 
Eq. (|4.88|) by parts with respect to X\,y\ and x 2 ,y 2 and using 

A X1 A yi 'y BA (x, y; x x y x ) = 5 (x - x x ) 5 (y - y x ) (4.89) 

we find that function 4>(x, y; 00) satisfies the equation 

= J d 2 zK(x, y;z) ((<f>(x,y;oo) - <p(x, z; 00) <p(z, y; 00)) (4.90) 

To obtain Eq. (|4.90|) from Eq. (|4.88|) we need to assume that 

J d 2 z K(x , y; z) A x A y 0(x, y; 00) « A x A y J d 2 z K(x, y; z) <f>(x, y; 00) (4.91) 
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Indeed, as we show below A x A y </>(x, y; oo) is very singular and behaves as A x A y <p(x, y; oo) oc 
5 (x — Xi) 5 (y — j/i). The function / d 2 z K(x, y; z) is less singular. Therefore, as far as the most 
singular part of solution is concerned Eq. (j4.91|) holds. 

The detailed discussion of the solution to Eq. f)4.88|) will appear in a separate paper jl?] and 
here we will only use the fact that 4>(x, y\ oo) = 1 satisfies this equation. 

Therefore, the cumulant generating functional for the asymptotic solution is equal to 

= « / d 2 xd 2 y A x A y( p(x,y;Y) j(x,y) (4.92) 

Using Eq. (|4.92J) and unitarity constraint of Eq. (|2.22j) . we can find the asymptotic behaviour for 
the scattering amplitude. One can see that the answer is 

N(x,y;x',y';oo) = 1 - (4.93) 

- exp (^-k 2 J d 2 xd 2 yd 2 x" d 2 y" A x A v <f)(x,y,x, y, oo)^y BA (x, y-,x" ,y")A x » Ay>,<f)(x',y'-,x" ,y" ,oa) 
Using the explicit form for •y BA as well as Eq. ()4.89|) . we can rewrite Eq. ()4.93|) in the form: 

N{x,y;x',y';oo) = 1 - (4.94) 

exp (~^^ J d 2 x" d 2 y" <j>(x, y; x" , y n ; oo)) A x » A y ^(x' , y'; x" , y" ; oo)J 
where is the solution of Eq. f|4.9U|) . 



4.1.2 Approaching the asymptotic solution 

As in the toy model we will search for a solution in the following form: 

$(y;[ 7i ]) = $(oo;[7 4 ]) + A$(r;[ 7l ]) (4.95) 

assuming Eq. ()4.86|) . and neglecting the contribution oc (A$) 2 . 
The resulting equation for A$ (Y; [7^]) has the form 

d A$(F;[ 7i ]) 



J d 2 x } d 2 yd 2 z d 2 x\ d 2 y x d 2 X2 d 2 y2 (4.96) 



dy 

{u(x,y) - u(x 1 ,y 1 )u(x 2 ,y 2 )) ( K(x,y;z)S (x - x x - x 2 ) 5 (y - y x - y 2 ) - 2 A x A y K(x, y; z) x 

1 BA (x, z; xi y x ) A X1 A yi (f)(x 1 , y x \ 00) >y BA (x, z; x 2 y 2 )) 9 ^} Y ' 

1 o-f(x 2} y 2 ) 

Eq. (|4.96|) is a Liouville-type equation which can easily be solved assuming that the functional 
A$ (Y; [ji]) = A$ ([7i(3^; yi)]) . Using this assumption we can re- write 

9A$(y;[7i]) f , 2 , 2 dA<S>{y;[n\) d 1 {y;x 1 y) 



J d 2 x d 2 y 



dy J * d>y(y;x,y) dy 
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and reduce Eq. (|4.96|) to the following form: 

dj(y;x,y) 



(4.97) 



dy 

d 2 z K(x, y; z) (u(y; x, y) - u(y; x, z) u(y; y,z)) -2nA x A y J d 2 zd 2 x 2 d 2 y 2 K(x, y; z) x 

7 BA (x, z; x 2 , m) <p(x, z; oo) ( u(y; x, y)) - uty x 2 , y%) u(3>; x, z)) 

where function u = 1 — 7 and 7 (y-, x, y) are determined by the initial condition at Y — 0, namely, 
at Y = we have only one dipole with coordinate x, y with 

A$Cy = 0;[ 7i ]) =ln(l- 7 (y = 0,x,y)) + k j(y — 0, x, y) A x A y <f>(x, y, 00) « (4.98) 

~ - 7(3> = 0, re, j/) + K>y(y = 0,^,2/) A x A y (p(x,y, 00) = II(a;,2/) 7^ = 0,x,j/) 

Using Eq. (J4.98j) we can re-write Eq. (J4.97j) in terms of A$ (y = 0; [7?]). Neglecting terms of 
order (A$) 2 we obtain the following equation 

= (4.99) 
y d 2 z (^K(x,y;z) - 2nA x A y J d 2 x 2 d 2 y 2 K(x, y; z) j BA (x, z; x 2 , y<i) (f>(y, z; 00) J x 

(2U(x,y) A ^'' X [ z) - *0>;x,y)) 
\ U{x,z) J 

Eq. ()4.99|) can be solved in log approximation. The first term in this equation is the familiar 
BFKL linear equation, therefore, we need to evaluate the second term in Eq. ()4.99j) . 

First, we assume that at high energies the nonlinear corrections set the new scale : the satu- 
ration momentum Q S {Y) PI El El and the integration which we have in the equation is really cut 
off at \x — z\ and /or \y — z\ smaller than 1/Q S |46| I48j. Actually, <f){y,z; 00) reaches a maximum 
value (f>(y, z\ 00) = 1 which is a solution to Eq. (|4.90p . We can rewrite K(x, y; z) in the form 

K(x,y;z) = A K ]n^ U ~ 



(y- K 



2 



Integrating by parts with respect to z, and differentiating j BA with respect to x only, we find that 
the second term has the form 

BA/-? ~ .-^ - s on/. .A^l^V) 



2 k A x A y J d 2 x 2 d 2 y 2 K(x, y; z) -f BA (x, z; x 2 , y 2 ) (p(y, z; 00) 2 U(x, y) 



U(x,z) 



2 K,A y ( d 2 x 2 d 2 y 2 \n(^^) A x A zl BA {x , z; x 2 , y 2 ) 2 H(x, y) X2 'f ) = (4.100) 

J \{y-z) J n(x 2 ,y 2 ) 

= 2 f d 2 z 1 A$(y;x,z) 
J (y- z) 1 
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where we put U(x,y) = U(x,z) to within the logarithmic accuracy. 

Collecting both terms and using Eq. ()4.89|) . we obtain the following equation for A$ 

d A$ (V-x v) r n 1 

fj- 211 = - J d 2 z-^—^(2A$(y;x,z) - A$(y;x,y)) (4.101) 

This is the BFKL equation. However, we solve this equation assuming that In [(a?— y) 2 Q 2 .)] S> 1. 
Considering \x — z\ <C \x — y\ we reduce Eq. (I4.1(jlj) to the following equation 

dA$(y-,x,y) f(*-y) 2 d 2 [x - z) 



d y Jp 2 (x — z) 

where A$ (y-, x, y) = A$ (y-, x, y) /{x - y) 2 . 

Actually, the typical cutoff p is of the order of p 2 = 1/Q 2 {Y) in the saturation region jUJJHB]- 
Eq. (|4.1U2j) can be re-written in the differential form 

* (4.103) 



dY d ln((x - y) 2 Q 2 (Y)) 
-A$(y-x,y) - ln((x - yf Q 2 S (Y)) 



. dA$0>;x,y) 



d \n((x - y) 2 Q 2 S (Y)) 



To find the solution to Eq. (j4.103|) we change the variable introducing a new variable 

HUES] 

z = ln((x - y) 2 Q 2 (Y)) = a s f^-Y + ln((f - y) 2 Q 2 ) (4.104) 
where j cr is determined by the following equation [TJ H5J HI] 

X(7cr) _ d Xilcr) 

1 ^y cr c? ^y cr . 

In terms of the new variable, Eq. ()4.103|) has the form 



(4.105) 



f^l = U {z) _ ,^M) ,4.106) 

and the solution of Eq. ()4.106j) is 

A$ = U{x,y)i{x,y) exp f- i-Z-^ ^ (4.107) 

To our surprise the asymptotic behaviour of A<3> turns out to be the same, as the asymptotic 
behaviour of the Balitsky-Kovchegov equation |4*Bj . 
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Using the unitarity constraint of Eq. (|2.22j) . and repeating all the calculations that led us to 
Eq. ()4.94|) we obtain, that the main corrections to the asymptotic behaviour given by Eq. ()4.94|) 
have the form: 

N{x,y;x',y',Y) = N (x, y; x', y', oo) (l - J d 2 x" d 2 y" U(x,y;x" ,y")Il(x',y';x" ,y") x 

( exp (" 5s S) {Y - rr ) + exp (- as S) (r)2 ))) (4 - 108) 

We assume in Eq. (|4. 108)1 that Y' and Y — Y' are much larger than hx(x — y) 2 Ql or \n(x' — y') 2 Ql 
. For large Y the minimal corrections occur at Y' = ~Y . Therefore 

R = ^i^i^j^jZj ~ jv* to y; ^, 2/'. oo) K exp / q a 1 ~ 7cr (y) 2 ^ (4109) 
N K (x,y;x',y',oa) \ S 8 x(Tcr)) / 

where we denote by k the fact that we consider the example with r(2 — > 3) = 0. 
It is interesting to compare this result with the solutions that have been found : 

• For the Balitsky - Kovchegov equation the ratio of Eq. (|4.1U9J1 is of the order of 

fl=Kexp (-" s 5S (i ' )2 ) 1 (4 - U0) 

• For Iancu-Mueller j2H] approach which can be justified only in the limited range of 
energies we have 

*^ exp (- as iS (y)2 ) ; (4111) 

Comparing Eq. ()4. 110)1 and Eq. (|4. 111)1 with Eq. (|4. 109)1 we see that the process of merging of two 
dipoles into one dipole crucially change the approach to the asymptotic behaviour . 



4.1.3 Main results 

The main results of our approach are given by Eq. (j4.108j) and Eq. (j4.109j) . These equations 
show that at high energies the scattering amplitude approaches the asymptotic solution which 
is a function of coordinates (N (x, y; x', y', oo), see Eq. (14.94)) ). This function is smaller than 1 
(N (x, y; x', y', oo) < 1) as it should be due to the unitarity constraints 

The corrections to the asymptotic behaviour turn out to be small (see Eq. (|4.109jl ). They are 
suppressed as exp(— CY 2 ), but the coefficient C is found to be 4 times less than for the solution 
of the Balitsky-Kovchegov equation. 

The toy model has proved to be very useful, and we showed that this model can serve as a 
guide, since we reproduce all the main features of the toy model solution, in the general solution 
to the problem. It should be stressed that the semi-classical approach based on a large parameter 
K ^> 1, gives the correct approximation to the problem. 
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4.2 General consideration 

4.2.1 Solution at Y -> oo for T(2 -»■ 1) = 

As in section 4.1.1 we search for the asymptotic solution, as the solution to the master equation 

(see Eq. (12.25)) ) with zero l.h.s. Assuming Eq. (|4.8fi(l and using Eq. for k we try to solve 
Eq. ()2.25|) looking for a solution of the form 

ihf r n ~i ( I d2xd2 V gpM/) l{x,y) \ (a-\-\o\ 

$(oc, M) = - k In f /rf2xrf2ye(X;j/) J (4-112) 



k In 



'^1 
.JV, 

where G is a function which is equal to 1 for |x| < i? and |y| < R and O = for \x\ > R and 
|y| > i? where i? is the largest scale in the problem, and 

^■(oo,[7i]) = f d 2 xd 2 y Q(x,y) j(x,y) and iV = J d 2 xd 2 y Q{x,y) (4.113) 
Substituting in Eq. (|2.25j) we obtain the following equation 

/© (x y) 
d 2 xd 2 y d 2 zT((x,y)^ (x, z) + (z, y)) ( 7(2;, y) - 7(2, z) 7(2, y)) — y 

k 2 J d 2 x x d 2 yi d 2 x 2 d 2 y 2 d 2 z T((xx, y x ) + (x 2 , y 2 ) -> [x 1: z) + (z, y 2 ) + (x 2 , yx)) x 

(7(^1, 1/2) 7(^2,2/1) ~ 7(gi» *) 7(^> 2/2) 7(^2, yi)) ^ Q (4.114) 

Since r((xi, yx) + (x 2 , 2/2) -> (zi, 2) + (2, y 2 ) + (x 2 , J/i)) = « T((xi, j/ 2 ) -> (xi, 2) + (2, y 2 )), these 
two terms cancel each other if we denote the variable of integration in the first term as (x x , yx) or 
as (xx,y 2 ) instead of (x,y). 

From Eq. (|4.112j) we obtain the generating functional at Y — > 00, which is equal 

Z(XM) = (S***V Q (*>vW*>v)y* (4.115) 

One can see that Eq. (|4. 115(1 cannot be normalized by the condition Z(Y, [uj = 1]) = Z(Y, p% = 
0]) = 1. Indeed, as for the case of the toy model (see Eq. (|3.6Uj) ) Z(Y, [%]) diverges at small 7's 
and the value of j 0ti should be fixed by including T(2 — ► 1) ^ 0. Including T(2 — > 1) we obtain 
the same normalization condition, but at 

li = 7o,i ~ «/« < 1 ■ 

The solution that satisfies this modified initial condition (Z(Y, [7^ = 70,1]) = 1) has the form 

7(v\ w ( I d 2 xd 2 yQ(x,y)'y( y x,y) \~ K 

Z{YM) = { %M J (4U6) 
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Using Eq. (|2.2U|) and Eq. 1)4.116)1 one calculates p n which are 

*(*,,»;... ») = nHfi^^ (4 . 117) 

We will find 70,* i n the next subsection. 

Using Eq. (|2.2()jl and the unitarity constraint of Eq. 1)2.22)1 we can find the asymptotic behaviour 
for the scattering amplitude. 

N((x,y;x',y';w) = (4.118) 
/ „ J d 2 xd 2 yd 2 xd 2 y Q(x,y) Q(x,y) 7 SA (x,y;x,y)\ 
~ 2 V ' ^ 5 ^ 2 (oo; [70,]) J 

where 2-P0 is the generalized hypergeometric function [55] , 

To obtain the answer we need to estimate the value of 70. 
4.2.2 7o(aj,y) 

As we have learned from the toy model, we need to compare the second term of Eq. 1)2.25)1 or 
Eq. ()2.11j) . with the third term in these equations using the asymptotic solution of Eq. 1)4.112)1 . 
7o(x,y) can be determined from the equation which equates the second term to the third one. 
Namely, we have 

d 2 x d 2 yy (x,y) T 2 ^i {{x,z) + (z,y) — > (x, y))d 2 z = J d 2 x 1 d 2 y 1 d 2 x 2 d 2 y 2 (4.119) 

7o(^i, 2/2) 70(^2, yx) r 2 ^ 3 ((x 1 , yx) + (x 2 , y 2 ) -> (2, Vi) + (ar 2 , z) + (^1, 2/2)) d 2 ^ 

Using the explicit forms of the vertices (see Eq. ()2.9j) and Eq. ()2.18)1 ) we can re- write Eq. ()4.119)) 
in a more convenient form, namely 

-J d 2 xd 2 y l0 {x,y) A x A y J d 2 z K{x,y;z) ^ BA {x, z; x lVl ) -y BA (z, y; x 2 y 2 ) = (4.120) 

= T^r 70(^1,2/2) lo(x 2 , yi) f (K(x 2 ,yx;z) + K(x l ,y 2 \z))d 2 z 

Taking A Xl A yi A X2 A y2 from both sides of Eq. ()4. 120)1 and integrating over Xi,x 2 ,yi,y 2 we can 
reduce this equation to a simpler form using the following properties of the Born amplitude 

A x A yl BA (x,z;x iyi ) = S^(x-xx)8^(y-yx) + 6®(iS - yx) 5^{y - xx) (4.121) 

Indeed, we reduce Eq. (14.120)1 to the form 

- J d 2 xd 2 y A x A y j (x,y) Jd 2 zK(x,y;z) = (4.122) 

= -= J d 2 x 1 d 2 y 1 d 2 x 2 d 2 y 2 A Xl A y2 -f (x 1 ,y 2 ) A X2 A yi (^j (x 2 ,yx) J d 2 z K (x 2 , yx] z) 
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Approximating 

A X2 A yi (j {x2,yi) J d 2 z K (x 2 ,yi;z)^j « A yi j (x 2 , yx) x (J d 2 z K (x 2 , yn z)) (4.123) 



we obtain that 



/ d 2 xd 2 y A x A ylo (x,y) = ? (4.124) 

J K 



The solution to Eq. (|4.124j) is 



lo(x,y) = -^— 2 - ln 2 ^. (4.125) 
16 n A k y 1 



One can check that Eq. ()4.123|) is valid for 70 (x, y) given by Eq. ()4.125|) at least as far as the most 
singular terms are concerned. 

Using Eq. ()4.125|) we can calculate the ratio 



J d 2 x d 2 y d 2 x d 2 y 0(x, y) 9(x, y) r ) BA (x, y; x, y) 
F 2 (oo; [7o,i]) 

First we calculate JF(oo; [7^] rewriting Eq. (J4.125j) in the form of contour integral 



(4.126) 



, v 1 1 f a + io ° dh ( x 2 \ h 

7o*,V = —— I — - (4.127) 



167T 2 2ni Ja-ioo h 3 \y 2 
Using Eq. (|4.127|) we can reduce Eq. (|4.113|) to the equation 

Hoc; ho,} = -t^tt- TF / dt(t h + t~ h ) = (4.128) 

k 16 2m Ja-ioo h 6 J Jo 

R 1 R 4 1 r a +^ dh 1 iifi 4 
k16Y 2m A-ioo 7^ 1 - h 2 k 32 
In Eq. ()4.128|) t is the new variable t = x 2 /y 2 . 

We use the same trick to calculate the numerator of Eq. (|4.126|) . We, finally, obtain 

/ d 2 x d 2 y d 2 x d 2 y Q(x,y) Q(x,y) ■y BA (x,y; x,y) 32k 



F 2 (oc; [7o,i]) 



l-2 



(4.129) 



Therefore, we have 



32 k n 

N(x,y;x',y';oo) = 2 F («,«;; -— ) < l-e" 32 * (4.130) 
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From Eq. (|4.13U|) we see that the amplitude turns out to be less that unity, and the unitarity 
limit cannot be reached. Therefore, we have gray disc scattering instead of the black disc one 
which we expect. 

4.2.3 Corrections to the asymptotic solution 

Substituting Eq. (|4.95jl into Eq. (|2.25jl with T(2 — > 1) = we obtain the equation for A$, namely, 



+ 2 1 — - f d 2 zd 2 x 1 d 2 y 1 d 2 x 2 d 2 y 2 K{x 1: y 2 ;z) x (4.131) 

./MOO, I'Ynfl ) J 

5A$(Y; [ 7i ]) 



F(oo, [7^) 

7(^1,2/2)7(^2,2/1) - 7(^1,2)7(^,2/2)7(^2,2/1)) 



^ 7(^1, 2/i) 



This Liouville-type equation can be solved assuming that A$(F; [ 7 j]) = A<&([ 7 j(y; a^,^)]) as 
we did in section 4.1.3. For function ^(Y; x, y) we obtain the equation 



d xd y K(x,y;z) (-f(y-,x,y) - 7(^5 x, z) u(y; z,y)) + 



dy 

+ 2 ^f/ d 2 zd 2 x 1 d 2 y 1 d 2 x 2 d 2 y 2 K(x 1} y 2 ',z) x (4.132) 

(70^1, 2/2) 7(^2, 2/1) - 7(3 ; ;xi,2)M(3 ; ;2,?/ 2 )7(y;x2, j/i) ) 

It is more convenient to switch to the funcion 7: m(3^;x,?/) = 1 — j(y;x,y). Indeed, even from 
the asymptotic solution we see that the typical values of u should be close to unity, namely, 
1 — w = 7 m 1/k. The relation between functions u and/or 7 is determined by the initial 
condition, namely, Z(y-, [itj) = u(x,y) at Y — 0. For the functional A<I>, this condition can be 
written 



A$(;y = 0; [ 7i ]) = ln(l- 7 (y = 0;x,y)) + k In 



7 d 2 xd 2 yQ(x,y) 7 (y = 0;x,y)' 
^(00; [ 7o ,i) 



~ - tCV = 0; rr, y) - win — = — (4.133) 

^(00; [ 7o ,ij) 

where A7(^; x, y) = 7(3^; 2/) — 7o(3^; a;, 2/)- The key idea of our approach to corrections at high 
energies, is to take A7 <C 70. In the derivation of Eq. ()4.133|) we have used this assumption. 

We can also write the relation for / d 2 xd 2 y A$(^ = 0; [7^]), namely, 

J d 2 xd 2 y A$CV = 0;[ 7i ]) = ^(00; [ 7o /J) - J d 2 x d 2 y Q(x, y) A j(y = 0; x, y) x (4.134) 

/ N \ 

1 - k 



F(oo; [ 7o , i ]) j 
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Using Eq. ([4.133)1 and Eq. ([4.13411 we obtain that 

1 (y = Q;x,y) = - A^(y = 0;x,y) + — " x (4.135) 

F(oc; [7o,iJ) + k N 



J d 2 xd 2 y A$(y = 0; x, y) + ^(oo; [ 7o , 4 ]) 



Introducing 7 (y-, x, y) = A(y-,x,y) + 70(00; a;, y) and considering A (y-, x, y) <C 70(00; x,y) 
we rewrite Eq. (14.132)1 as a linear equation with respect to A (y-, x, y), namely 



dy 



J d 2 zK(x,y;z) (Aj(y;x,z) + A 7 (y-z,y) - A 7 (y-x,y)) (4.136) 



Eq. (|4. 136)1 is the BFKL equation, but with extra sign minus in front. Actually, this is the 
same equation as Eq. ()4. 101)1 which has been solved (see Eq. ([4.107)1 ) . The solution can be written 
in the form: 

A^y-x^y) = A 7 (y = 0;x,y) exp (- i "J^ z 2 ^j (4.137) 

where z is defined in Eq. ([4.104)1 . 

The initial condition of Eq. ([4.135)1 determines the relation between function 7 and the func- 
tional A$. It is even more convenient to use the initial conditions in the form of Eq. ([4.135)1 or 
even of Eq. (ITOHj) 

Using the unitarity constraint of Eq. (|2.22|) . and following the same line of calculation as in 
the derivation of Eq. ([4.118)1 we obtain: 

N(x,y-x\y\Y) = N (x, y; x\ y\ 00) (1- (l + «)x 

K- as S) (Y - rr ) + exp (- as S) (r)2 ))) (4 - 138) 

We assume that Y' and Y — Y' are much larger than ln(x — y) 2 Q 2 . Once more the minimal 
corrections occur at y' = |y and we have 

N (x,y;x',y',Y) - N R (x, y; x', y', 00) / _ 1 - 7 cr 2 \ . . 

R * = N k(x , r , xW ^) K exp I " as s^T) (Y) ) (4139) 

By index k we denote the fact that we consider the case with r(2 — > 1) = 0. N% (x, y; x', y', 00) is 
given by Eq. (14.130)1 . 

In spite of the fact that the asymptotic solutions look quite different, both examples approach 
their asymptotic behaviour at the same rate. 



31 



4.2.4 General property of the solution 

Eq. ()4.139|) shows that the general solution has the following attractive features 



• The asymptotic solution is not unity, but a function of coordinates which is smaller than 
1. This means that the high energy scattering corresponds rather to the gray disc regime 
and not to the black disc one which was expected; 

• The approach to the asymptotic solution is a very fast decreasing function of energy which 
is proportional to e~ CY . The coefficient C in 4 times smaller that for the solution [IS] of 
the Balitsky-Kovchegov equation and in 2 times smaller than for the Iancu-Mueller solution 

B3|. ' 

• The entire dynamics can be found, neglecting the process of the merging of two dipoles 
into one dipole. This process we need only to determine the value of 7 , but it does not 
affect any qualitative properties of the asymptotic behaviour of the scattering amplitude at 
high energy. 

In general we found the asymptotic solution to the evolution equation (see Eq. (|2.11j) ) and 
we prove that this solution is self consistent. It means, that corrections to this solution decrease 
with energy. However, we did not prove that our solution satisfies the correct initial conditions for 
the generating functional. Therefore, in principle, we could overlook a solution which is diffrent 
from ours but satisfies the initial condition while ous does not satisfy them. Such scenario looks 
unlikely to us because of our experience with the toy model and because of the exact solution 
to Eq. ()3.32j) . These examples show that we have a correct procedure for finding the asymptotic 
solution. 

5 Conclusions 

In this paper we found for the first time, the high energy amplitude in QCD taking into account the 
Pomeron loops. We found that the processes of merging between two dipoles as well as processes of 
transition of two dipoles into three dipoles crucially change the high energy asymptotic behaviour 
of the scattering amplitude. 

Using the fact that in QCD, T(l -> 2)/T(2 -»• 1) > 1 as well as T(l -> 2)/T(2 -> 3) > 1 we 
developed a semi-classical method for searching the solution. The main results of our approach are 
presented in Eq. (jOi|) .En. (|4.1t)9|) . Eq. (I4.13t)j) and Eq. (|4.139|) . They show several unexpected 
features of high scattering in QCD: (i) the asymptotic amplitude is a function of dipole sizes and, 
therefore, the scattering amplitude describes the gray disc structure at high energy, instead of black 
disc regime which was expected; (ii) the solution approaches the asymptotic in the same way as 
the solution to the Balitsky-Kovchegov equation (oc exp(— CY 2 ) ) but coefficient C in four times 
smaller than for the Balitsky-Kovchegov equation; (iii) the process of the merging of two dipoles 
into one only influences the initial condition for the generating functional. Indeed, the new initial 



32 



conditions is Z(Y; [ui = 1 — 7o,i]) = 1 instead of Z(Y; [it, = 1]) = 1. The value of 7 is determined 
by the process of the merging of two dipoles into one, and it turns out to be of the order of 
k/n <C 1. Putting this new initial condition we can use the Balitsky-JIMWLK approach [T4*l ITS] 
for description of the high energy asymptotic behaviour of the scattering amplitude without any 
modifications recently suggested [TH] . 

One of the experimental consequences that stems from our solution, is the fact that we do not 
expect geometrical scaling behaviour of the amplitude in the saturation region. We also do not 
expect the black disc behaviour, which in the past was considered as the most plausible for the 
high energy asymptotic behaviour. 

In general, the generating functional approach has many advantages, in particular, it is very 
simple with a clear control of physics in each step of evolution. However, it has its own shortcoming: 
inside of this approach we cannot restrict ourselves by the finite nuimber of possible transitions 
5 . For example, why we included T(2 — > 3) and did not include T(3 — > 2)? Formal argument 
is that we included all vertices of the order of l/iVf and neglected all vertices of the order of 
1/iVg. However, without T(3 — > 2) we do not have symmetric approach and, strictly speaking, 
we cannot guarantee the t-channel unitarity. Since we use t-channel unitarity in our estimates 
for the scattering amplitude, we obtain the answer which respects the t-channel unitarity, but, 
in the spirit of Iancu-Mueller factorization (22], this answer is correct only for limited range of 
energy. The work is in progress for taking into account T(3 — > 2), however, the final answer to the 
question which vertices we need to include, perhaps, lies in more general theoretical approaches 
(see Refs, [HI HH HQ 122 Eli ) • 
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